Stennis Space Center 



Assessing MODIS-based products and 
techniques for detecting gypsy moth 

defoliation 

Presented by Joseph P. Spruce 
Science Systems and Applications, Inc. (SSAI) 

John C. Stennis Space Center 

Co-Authors: 

William Hargrove, Robert E. Ryan, James Smoot, Don Prados, Rodney McKellip, 

Steven A. Sader, Jerry Gasser, George May 

2008 US-IALE Meeting 
Madison, Wisconsin - April 9, 2008 


RELEASED - Printed documents may be obsolete; validate prior to use. 


Stennis Space Center 


Project Background 

• The project addresses: What is the potential of MODIS data 
for monitoring historic gypsy moth defoliation? 

• This project regards a NASA/USDA Forest Service (USFS) 
partnership 

- NASA is helping USFS to implement satellite data products 
into its emerging Forest Threat Early Warning System 

- The latter system is being developed by the USFS Eastern 
and Western Forest Threat Assessment Centers 

- The USFS Forest Threat Centers want to use MODIS time 
series data for regional monitoring of forest damage (e.g., 
defoliation), preferably in near real time 
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Initial Study Area - 
MickAppalachian Highlands 


Study Area Outlined in Yellow Below (Total -15.5 Million Acres) 
Gypsy Moth Defoliation Occurred on Multiple Years During MODIS Era 




Note: This Area Encompasses Several Landsat Scenes 
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Project O bjecti ves 



1 . Assess 250-meter MODIS time series data for mapping 
historic gypsy moth defoliation 

- Assess accuracy of detection products compared to reference data 

• Initial focus on 2001 due to known extensive defoliation and available 
reference data 

- 2001 Landsat and ASTER imagery 

- 2001 USFS defoliation sketch maps 


2. Assess simulated VIIRS time series data for mapping 
historic gypsy moth defoliation 

- VIIRS will be the follow-on to MODIS 

- VIIRS will have 400-meter resolution 

- VIIRS data is being simulated mainly from MODIS imagery 
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Reducing Noise in the 
MODIS and VIIRS Time Series Data 



• Data from each MODIS-based product was independently 
preprocessed to reduce inclusion of clouds and other low- 
quality data 

- MODIS MOD02 (planetary reflectance data) and MODI 3 
(atmospherically corrected NDVI) 

• The Time Series Product Tool (TSPT) software was used 
to derive “cleaned” MODIS and VIIRS time series data 

- TSPT outputs vegetation index products (e.g., NDVI) 
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Computing Maximum NDVI Composites 
for theHefoliation Time Frame 



• Computed a maximum NDVI composite for the gypsy moth 
defoliation period of each year (June 1 0— July 27) 

- Time series includes defoliated and non-defoliated years 

• Stacked maximum NDVI during defoliation image for 2001 
with maximum NDVI image during the defoliation period 
over the whole 2000-2006 time series 

• Used data stack to compute defoliation detection products 



MODIS Maximum NDVI Defoliation RGB - Based 

Defoliation Time Frame (DTF) Maximum NDVI DTF Maximum NDVI DTF on Maximum NDVI DTF 

for 2000-2006 (Red) 2001 (Green) 2001 (Blue) for All Years vs. 2001 
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Views of Gypsy Moth Defoliation on 
MODIS versu s Landsat NDVI D ata 
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MODIS - Maximum NDVI DTF All Years in Red; Same for 2001 in Blue and Green 

DTF - Defoliation Time Frame (June 1 0— July 27) 

Landsat - 6/1 0/2000 NDVI Loaded in Red; 7/15/2001 NDVI Loaded into Blue and Green 



• Both RGB images show defoliation from 2001 outbreak in red tones 

• MODIS RGB is cloud free due to temporal processing of daily data 
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Method for MODIS Image Classification 
of Gypsy M oth D efoliation 



• Processed 2001 MODIS and VIIRS data into defoliation 
maps of 2 classes: defoliated versus other 

- Employed unsupervised classification techniques 

- Example results from MODI 3 250 m, MOD02 250 m, and simulated 
VIIRS 400 m data 

• Applied post-classification “filtering” technique to reduce 
commission errors from patches smaller than 1x1 km 


Defoliation RGB - Based 20-Class Unsupervised Application of 1x1 km Final Result 

on Maximum NDVI DTF for Classification Clump Aggregation Overlain onto 

All Years vs. 2001 Defoliated Forest in Red Filter NLCD Forest Mask 



I^^OTa/^rorauti^amf^paS^tcfminls^ration 3 ^ obsolete, validate prior to us ^ ssess j n g MODIS-based products and techniques for detecting gypsy moth defoliation 8 






2001 Defoliation Classifications from MODI 3, 
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MODG2, and Simulated VIIR S NDVI Pr oducts 

er 7 
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Method for Accuracy Assessment 
of_Def o I i a t i orLDetecti o n_Ero d u ct s 



• Drew stratified random sample locations from best apparent 
classification (MOD02 250-meter result) 

- Drew samples for defoliated forest versus other classes 

• An image analyst interpreted each sample location Landsat 
or ASTER as being defoliated or other 

• Interpretation results were then compared to each test 
classification 

- Examples - MOD02, simulated VIIRS, and MODI 3 products 

• Final results were summarized for defoliated forest versus 
“other” class 
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Relative Accuracy of Example 2001 
Defoliation Classification Products 



2001 Classification 
Product 

Defoliated Forest 

Other 

Overall 

PA 

UA 

Kappa 

PA 

UA 

Kappa 

OA 

OK 

MOD02 NDVI 250 m 

91 % 

( 52 / 57 ) 

78 % 

( 52 / 67 ) 

0.67 

87 % 

( 101 / 116 ) 

95 % 

( 101 / 106 ) 

0.86 

88 % 

( 153 / 173 ) 

0.75 

VIIRS NDVI 400 m 
(Simulated from 
MOD02) 

86 % 

( 49 / 57 ) 

78 % 

( 49 / 63 ) 

0.67 

88 % 

( 102 / 116 ) 

93 % 

( 102 / 110 ) 

0.78 

87 % 

( 151 / 173 ) 

0.72 

MODI 3 NDVI 250 m 

44 % 

( 25 / 57 ) 

86 % 

( 25 / 29 ) 

0.79 

97 % 

( 112 / 116 ) 

78 % 

( 112 / 144 ) 

0.33 

79 % 

( 137 / 173 ) 

0.46 


Note: PA = % Producer’s Agreement (# correct/total), UA = % User’s Agreement (# correct/total), 
Kappa = Kappa Statistic, OA = % Overall Agreement (# correct/total), and OK = Overall Kappa. 
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Conclusions for Example 
2001 def oliation Mapping Products 



• MODIS and simulated VIIRS time series data produced 
effective regional defoliation maps for 2001 

- Temporal processing techniques and pest phenology 
knowledge aided the application 

• MOD02 daily products yielded the best results 

- MOD02 250 m and simulated 400 m VIIRS NDVI products 
yielded similar measures of accuracy 

- MODI 3 NDVI defoliation maps showed the lower overall 
accuracy, in part from omission of defoliation areas 
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Gypsy Moth Defoliation of 2007 


Preliminary Results Using Available MODI 3 NDVI Data 
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View of Entire MODI 3 Mosaic 
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Used in 2007 Case Study 


Maximum NDVI 2007 DTF in Red; Maximum NDVI DTF All Years in Blue and Green 
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View of 2007 Gypsy Moth Defoliation 
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From MODI 3 Data 


Maximum NDVI 2007 DTF in Red; Maximum NDVI DTF All Years in Blue and Green 


Total Land 
Area Shown 
~ 44 Million 
Acres 
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View of 2007 Gypsy Moth 
Defoliation Map from MOD 13 Data 
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Defoliation Map in Foreground is Overlain onto USFS 250 m Forest/Non-Forest Map 
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View of 2007 Gypsy Moth Defoliation 
From USFS Sketch Mao 


Sketch Map in Foreground is Overlain onto USFS 250 m Forest/Non-Forest Map 
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Relevance of Pr oject 



• The project showed potential of MODIS and VI IRS time 
series data for contributing defoliation detection products to 
the USFS forest threat early warning system 

• This study yielded the first satellite-based wall-to-wall 2001 
gypsy moth defoliation map for the study area 

• Initial results led to follow-on work to map 2007 gypsy moth 
defoliation over the eastern United States (in progress) 

• MODIS-based defoliation maps offer promise for aiding 
aerial sketch maps either in planning surveys and/or 
adjusting acreage estimates of annual defoliation 

• More work still needs to be done to assess potential of 
technology for “now casts” of defoliation 
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MODIS^Related Terminolog y 



• MODIS - Moderate Resolution Imaging Spectroradiometer 

- 2 MODIS instruments in space (Aqua and Terra satellites) 

• Each sensor collects 1 image per location each day 

- 250-meter resolution for NIR (near infrared) and red bands 


• NDVI - Normalized Difference Vegetation Index 

- NDVI = (NIR - red) / (NIR + red) 

- Can be computed from multiple MODIS products 

• MOD02 - Daily MODIS Radiance Data 

• MODI 3 - MODIS 16 Day Composite Vegetation Indices 

• VIIRS - Visible/Infrared Imager/Radiometer Suite 

- future follow-on to MODIS 
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Gypsy Moth Defoliation Maps from MODIS 
and Simulated VIIRS Da ta 


0 


MODIS (MOD02) and VIIRS Products Produced Similar Results 
Both Tended to Map Same Outbreak Areas 



Red - Defoliation 

Green - Forest from NLCD 

Tan - Non Forest from NLCD 
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Aerial View 

of Gypsy Mot h Defo liation 
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